Synovial fibroblast gene expression is associated with sensory nerve growth and pain in rheumatoid arthritis.

It has been presumed that rheumatoid arthritis (RA) joint pain is related to inflammation in the synovium; however, recent studies reveal that pain scores in patients do not correlate with synovial inflammation. We developed a machine-learning approach (graph-based gene expression module identification or GbGMI) to identify an 815-gene expression module associated with pain in synovial biopsy samples from patients with established RA who had limited synovial inflammation at arthroplasty. We then validated this finding in an independent cohort of synovial biopsy samples from patients who had early untreated RA with little inflammation. Single-cell RNA sequencing analyses indicated that most of these 815 genes were most robustly expressed by lining layer synovial fibroblasts. Receptor-ligand interaction analysis predicted cross-talk between human lining layer fibroblasts and human dorsal root ganglion neurons expressing calcitonin gene-related peptide (CGRP+). Both RA synovial fibroblast culture supernatant and netrin-4, which is abundantly expressed by lining fibroblasts and was within the GbGMI-identified pain-associated gene module, increased the branching of pain-sensitive murine CGRP+ dorsal root ganglion neurons in vitro. Imaging of solvent-cleared synovial tissue with little inflammation from humans with RA revealed CGRP+ pain-sensing neurons encasing blood vessels growing into synovial hypertrophic papilla. Together, these findings support a model whereby synovial lining fibroblasts express genes associated with pain that enhance the growth of pain-sensing neurons into regions of synovial hypertrophy in RA.

Machine Learning Reveals Synovial Fibroblast Genes Associated with Pain Affect Sensory Nerve Growth in Rheumatoid Arthritis.

It has been presumed that rheumatoid arthritis (RA) joint pain is related to inflammation in the synovium; however, recent studies reveal that pain scores in patients do not correlate with synovial inflammation. We identified a module of 815 genes associated with pain, using a novel machine learning approach, Graph-based Gene expression Module Identification (GbGMI), in samples from patients with longstanding RA, but limited synovial inflammation at arthroplasty, and validated this finding in an independent cohort of synovial biopsy samples from early, untreated RA patients. Single-cell RNA-seq analyses indicated these genes were most robustly expressed by lining layer fibroblasts and receptor-ligand interaction analysis predicted robust lining layer fibroblast crosstalk with pain sensitive CGRP+ dorsal root ganglion sensory neurons. Netrin-4, which is abundantly expressed by lining fibroblasts and associated with pain, significantly increased the branching of pain-sensitive CGRP+ neurons in vitro . We conclude GbGMI is a useful method for identifying a module of genes that associate with a clinical feature of interest. Using this approach, we find that Netrin-4 is produced by synovial fibroblasts in the absence of inflammation and can enhance the outgrowth of CGRP+ pain sensitive nerve fibers. One Sentence Summary: Machine Learning reveals synovial fibroblast genes related to pain affect sensory nerve growth in Rheumatoid Arthritis addresses unmet clinical need.

Computer Vision Analysis of Rheumatoid Arthritis Synovium Reveals Lymphocytic Inflammation Is Associated With Immunoglobulin Skewing in Blood.

OBJECTIVE: We sought to develop computer vision methods to quantify aggregates of cells in synovial tissue and compare these with clinical and gene expression parameters. METHODS: We assembled a computer vision pipeline to quantify five features encompassing synovial cell density and aggregates and compared these with pathologist scores, disease classification, autoantibody status, and RNA expression in a cohort of 156 patients with rheumatoid arthritis (RA) and 149 patients with osteoarthritis (OA). RESULTS: All five features were associated with pathologist scores of synovial lymphocytic inflammation (P < 0.0001). Three features that related to the cells per unit of tissue were significantly increased in patients with both seronegative and seropositive RA compared with those with OA; on the other hand, aggregate features (number and diameter) were significantly increased in seropositive, but not seronegative, RA compared with OA. Aggregate diameter was associated with the gene expression of immunoglobulin heavy-chain genes in the synovial tissue. Compared with blood, synovial immunoglobulin isotypes were skewed from IGHM and IGHD to IGHG3 and IGHG1. Further, patients with RA with high levels of lymphocytic infiltrates in the synovium demonstrated parallel skewing in their blood with a relative decrease in IGHGM (P < 0.002) and IGHD (P < 0.03) and an increase in class-switched immunoglobulin genes IGHG3 (P < 0.03) and IGHG1 (P < 0.002). CONCLUSION: High-resolution automated identification and quantification of synovial immune cell aggregates uncovered skewing in the synovium from naïve IGHD and IGHM to memory IGHG3 and IGHG1 and revealed that this process is reflected in the blood of patients with high inflammatory synovium.